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Abstract 

We present a systematic theory of dissipation in finite Fermi sys- 
tems. This theory is based on the application of periodic orbit theory 
to linear response of many-body systems. We concentrate only on the 
mesoscopic aspect of the phenomena wherein a many-body system can 
be reduced to a single-body in an effective mean-field. We obtain semi- 
classical periodic-orbit corrections on top of the two-time correlation 
function for the rate of energy dissipation. We show that this energy 
dissipation is irreversible on an observational time- scale. To do so, 
we derive a generalised Smoluchowski equation for the energy distri- 
bution in the quantal domain. Employing the Weyl-Wigner expansion, 
we also obtain the equation governing the evolution of energy distribu- 
tion in the combination of semiclassical and adiabatic approximations. 
Further, we show how the periodic orbit corrections are related to ge- 
ometric phase acquired by single- particle wavefunction as it evolves 
with the slow, time-varying mean-field. We present our results for the 
important case of mixed dynamics also. We obtain random-matrix re- 
sults for response functions. We incorpoate chaos in the underlying 
classical system by writing an ansatz for a generic wavefunction which 
leads to various central results of equilibrium statistical mechanics. 
This new formalism is extended here to include dissipation. Finally, we 
present an expression for the viscosity tensor encountered in nuclear 
fission in terms of periodic orbits of single particle in an adiabatically 
deforming nucleus. 

PACS Nos. : 21.60.-n, 24.60.Lz, 05.45.Mt 

Keywords : Chaos and dissipation, Quantum transport, Periodic orbit the- 
ory, Quantum diffusion and irreversibility 



*srjain@apsara.barc.ernet.in 



1 



1. Introduction 

Nuclei and metallic clusters support a number of collective excitations 
which are well-known and well-studied [2|. Damping of these collective 
excitations is not well understood. Dissipation related to chaotic motion of 
a single particle in an effective mean-field has been proposed and pursued 
for last two decades [||, We are interested in the mesoscopic aspect of this 
phenomenon where the mean free path of a single particle is of the order 
of the size of the system. For instance, the frequency at which an isovector 
giant dipole resonance is excited in a nucleus is such that the corresponding 
wavelength is larger than the nuclear size. In the standard parlance of nuclear 
physics, this is referred to as one-body dissipation. In this paper, we present 
a systematic semiclassical theory of one-body dissipation. Since it is well- 
known that several spectral properties of chaotic quantum systems, nuclei 
etc. are well represented by random matrix theory, we present results where 
random matrix theory is conjointly taken with linear response theory. 

This work is inspired by some recent developments where some of the 
most important results of equilibrium statistical mechanics have been de- 
rived, starting from chaos in the single-particle motion in gases. Due to 
presence of chaos, the eigenfunctions of the corresponding quantum system 
become irregular, where irregularity is quantified by some correlations. For 
instance, spatial correlation function of an irregular eigenfunction is given 
by a Bessel function || ^] . The eigenfunctions can be modelled by a random 
superposition of plane waves. Remarkably, from this ansatz, without thermal 
reservoir, all the momentum distributions (viz., Maxwell-Boltzmann, Fermi- 
Dirac, Bose-Einstein) have been re-derived by Srednicki ||. This approach 
plays a very important role in two dimensions where some of the standard 
methods of statistical mechanics fail. For the two-dimensional case, the mo- 
mentum distribution was found J5| for an anyon gas upto o(h 2 ) . It was 
also shown that the ansatz is consistent with random matrix theory. 
Employing a random matrix argument, the second law of thermodynamics 
could also be arrived at by Jain and Alonso |?j . In this work, we attempt to 
extend this framework to include one-body dissipation. 

Let us recall that for a nucleus participating in a reaction, the time-scale 
is ~ 10~ 20 s or even shorter. It has been qualitatively argued || that in 
this situation, a nucleus is an isolated system. Thus, we must arrive at the 
notion of dissipation for an isolated quantum system since we have already 
succeeded in arguing for thermalisation of isolated quantum many-body sys- 
tems ||. 

A formal theory of one-body dissipation was developed in classic works 
where linear response techniques were used ]Tl| , [l0[ ] . Our considerations also 
rest on linear response theory. Beginning from linear response theory, we 
present a semiclassical expression for response function. For sake of com- 
parison, we also present random matrix expression. From this, we obtain 
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the rate of energy dissipation. Clearly, if there is energy dissipation, the 
energy distribution must evolve irreversibly in the quantal problem. This 
becomes an important issue as there can be no irreversibility or dissipation 
in an isolated quantal system that supports only pure point spectrum |l^| . 
Employing multiple time-scale analysis, we show that for sufficiently long 
times, there is irreversibility if the spectrum is coarse-grained (e.g., by a 
random matrix ensemble average). Beginning from the von Neumann equa- 
tion for the density operator, we show that the energy distribution obeys 
a generalised Smoluchowski equation. When semiclassical and adiabatic ap- 
proximations are put together, it is a highly singular situation. We perform 
a double expansion, one in the small adiabaticity parameter, e, and the other 
in h. We then obtain the equation for energy distribution valid upto times 
of order e _1 , and upto order fi 2 . With this derivation, we have thus found 
the solution of the following general problem. If a quantal system is evolving 
adiabatically, the energy levels of the system will evolve in a complicated 
manner fl3|| . As the levels come much closer than the mean level spacing at 
a particular epoch of evolution, the non-adiabatic transitions become more 
probable. This can be understood by employing degenerate perturbation 
theory. Thus, if we distribute some particles initially (at time, t = 0) in 
accordance with certain law (e.g., the Boltzmann law), the final distribu- 
tion will be very different. It has been thought that the evolution of the 
energy distribution must be diffusive. Although the guess is qualitatively 
correct, the equation is not the Smoluchowski equation but another equa- 
tion which we have called as generalised Smoluchowski equation [14]. This 
general problem is encountered in many fields of physics and chemistry. 

Subsequently, we present relations connecting single-particle properties 
and bulk properties. We show that the geometric phase acquired by a single- 
particle wavefunction due to adiabatically vibrating mean-field is related to 
the response function. Although this result has been obtained earlier 15], 
we briefly discuss here to put it in proper perspective. 

All the above results assume that single-particle dynamics is chaotic. In 
realistic systems, all the symmetries are not broken. Thus the dynamics is 
usually mixed, i.e., there are stable islands embedded in stochastic sea. We 
obtain the results for this case which are based on trace formula for the case 
of partially broken symmetry obtained by Creagh ifTq] . 

The response function should depend on the shape of the cavity or the 
mean- field potential. In a relatively recent attempt, shape dependence was 
incorporated by multiplying by a factor, 71, the average fraction of trajecto- 
ries which are chaotic when a uniform sampling is done [17]. This is based 
on some heuristic arguments. It was shown that this chaos weighted wall 
formula is in a better agreement with experimental data [18]. This observa- 
tion and knowing that the spectral properties of non-integrable systems can 
be understood using semiclassical arguments based on the trace formulae 
[19 1, we believe that the theory developed here is logical and natural. Since 
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it is known that the spectral properties of chaotic quantum systems can be 
understood in terms of semiclassical trace formulae, our approach is more 
general and system-specific. We believe that the development presented here 
provides a better understanding of the interplay between single-particle mo- 
tion and collective motions in a many-body quantum system. 



2. One-body dissipation and wall formula 



To understand damping in low-energy heavy-ion collisions, the concept 
of one-body dissipation was introduced [|j where dissipation results from 
the interaction of nucleons with the time- varying mean- field. This involves 
a reduction of many-body system to a single-particle in an effective mean- 
field. Although two-body collisions are also important, due to Pauli exclusion 
principle, the available phase space for such collisions is much lesser and the 
one-body mechanism is expected to dominate in the mesoscopic regime. 

Assuming that the total one-body Hamiltonian is H(r, p; t) where (r, p) 
are phase space coordinates of a single particle. Splitting H into a time- 
independent and a time-dependent part, respectively Hq(t, p) and Hi(r, p; t), 
under the assumptions of linear response theory, the rate of energy dissipa- 
tion is given by fllTf , 

E = -J jr °di / fl-i(Bo(r,p;0,Po(r,p;t / );t) 

Fi(r,p;i)||(r,p), (1) 

where /o is the single-particle phase space distribution and an overdot rep- 
resents a time-derivative. The important quantity in (1) is the two-time 
correlation function in the integrand. 

The rate of energy dissipation for the case when the dynamics of a single- 
particle is chaotic is given by the wall formula jij : 

E W£L \\ = pv u 2 (a)d 2 a, (2) 



where u(a) is the normal component of the surface velocity at the point a 
on the surface, p and v are nuclear mass density and average nucleon speed 
inside the nucleus. 

The above expressions hold when the dynamics of a single particle is 
fully chaotic. However, since some of the symmetries are not broken, the 
dynamics is expected to be mixed rather than chaotic. For this case, a scaled 
wall formula [fijl], 

E = /Z-Evaii, (3) 

has been proposed. The factor, /Z, is the average fraction of the trajectories 
which are chaotic when a uniform sampling of the surface is done. The factor 
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is fixed from the parameter appearing in the Berry-Robnik distribution |2(J , 
which gives the nearest- neighbour level spacing distribution for systems 
with mixed phase space. This way of incorporating the role of dynamics 
is empirical. In our framework, we take care of the dynamical features ex- 
plicitly, hence any role of dynamics enters explicitly. Importantly, since we 
employ the trace formula, there is harmony between the dynamics and the 
spectral statistics. 

Let us now comment on the relation between chaos and dissipation. From 
a number of studies, classical and quantal f22|| , it is clear that if the dy- 
namics is chaotic, there is dissipation. If the particle and wall motions are 
considered self-consistently, then it has been numerically seen that the dy- 



namics is chaotic [23|. In this case, the energy is gained irreversibly, leading 
to dissipation. If the dynamics is regular or integrable, the energy gained by 
a particle from the wall will eventually be fed back to the wall, resulting in 
no dissipation. 

The one-body mechanism is dominant for low-energy nucleus-nucleus col- 
lisions or fission, i.e., in general, to slow collective processes. It is known that 
to simultaneously reproduce both the fission probability and the pre-scission 
neutron multiplicity, one needs a shape-dependent friction 24, 25]. These 



works came up as an attempt to combine dynamics and thermodynamics 
in a consistent way. The shape of the nucleus is reflected in the level den- 
sity and this is used, in turn, to construct entropy. At this point, whereas 
these works in nuclear physics accept a stochastic description in working 
with Langevin equation, we would persevere with the dynamical considera- 
tions. We believe that friction or viscosity, if any, in these quantum systems 
will have quantal signatures. For instance, the response function of a many- 
body Fermi system is related to geometric phase acquired by a single-particle 
eigenfunction as the system deforms [|15[| . These considerations need a sys- 
tematic theory where shape plays an important role. To do so, we believe 
that it is important to employ the Gutzwiller [26] trace formula and evolve 
a framework. 

3. Quantum chaos and equilibrium statistical mechanics 

The connection between chaos and equilibration or thermalization of 
many-body systems has been always believed in. The tenets of equilibrium 
statistical mechanics rest on the assumption of "molecular chaos" - an idea 
that goes back to Boltzmann [27]. With the advancement of our under- 



standing of the theory of dynamical systems and its relation with statistical 
mechanics [28|, it has become important to extend this understanding to 



quantal many-body systems. Only recently, such a connection has been es- 
tablished. Since this connection and the present approach is an underlying 
theme of this paper, we briefly present the status of our understanding for 
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want of a better perspective. 

Let us mention the numerical results on eigenfunctions of systems whose 
classical counterparts are chaotic. 

(Nl) The amplitude distribution of eigenfunctions are found to agree with a 



Gaussian distribution - we call these states as the generic ones [29]. Of course, 
we are not referring to eigenfunctions corresponding to the ground state 
which is special. Nevertheless, this observation is more general as the Gaus- 
sian distribution remains good even for pseudointegrable (non-integrable 
and non-chaotic) systems like a 7r/3-rhombus billiard [[?(], |3l| where the 
Kolmogorov-Sinai entropy is zero. 

(N2) The spatial correlation function, 

C a {s) = I J dxrfc (r - |) ^ (r + |) , (4) 

corresponding to an eigenfunction ip a ( r ) is found to agree with a (cylindrical) 
Bessel function [^], where V is volume. For two-dimensional billiards, it has 

been found to agree with Jo(ks), k = ^2mE/h 2 p9[j . 

(N3) The nodal lines are quite complicated, for pseudointegrable as well as 
chaotic billiards ||, ||. 



As far as eigenvalues or the energy levels are concerned, the measures 
quantifying spectral statistics agree with the results known from random ma- 
trix theory |[32|]. Even for pseudointegrable billiards, which are non-integrable, 
the nearest-neighbour level spacing distribution can be explained in terms 
of new ensembles (33] . 

The numerical results (N1,N2) can be explained if we assume that the 
complicated eigenfunctions are written as a random superposition of plane 
waves. This has been advocated by Berry || and written explicitly by Sred- 
nicki ||. For an iV-particle system, the canonical pair of coordinates are 
(X, P) where X = (xi, X2, xjv) and P = (pi, P2, Pjv)- The energy 
eigenfunction ?/> a (X) is thus written as 

^ a (X) = N a J d dN PA a (P)5(P 2 - 2mE a )e* x - p (5) 

for a d-dimensional system, where N a is the normalization constant and the 
amplitudes satisfy the two-point correlation function, 

xdN(-p _ p/\ 

(K(P)A,(P')) ME = S a , §( ^_ p/2 > . (6) 

This correlator has been shown to be consistent with random matrix the- 
ory, hence the appearance of a statistical average on the left side of the 
equation (ME stands for a matrix-ensemble average) ||. Beginning from 
this ansatz for the eigenfunction, the momentum distributions for ideal 
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gases have been shown to concur with the well-known results (like Maxwell- 
Boltzmann, Fermi-Dirac, Bose-Einstein, or fractional statistics) by Srednicki 
||, and, Jain and Alonso Even the first few virial coefficients have been 
calculated. The Wigner distribution corresponding to an eigenfunction is 
microcanonical. Thus, generically, corresponding to an eigenfunction, there 
is entire energy surface in phase space. 

An important question - in this theory, what is temperature ? Tempera- 
ture quantifies heat, and heat is defined as a mode of motion. Even if there is 
no reservoir, one can still consistently use the kinetic theory, and write for a 
system with / freedoms at energy E a , a temperature, T a , with E a = \fkT ai 
where k is the Boltzmann constant. As the average over the entire system 
involves an averaging over the level density, the average temperature, T will 
automatically get defined. This average temperature is the one we measure 
as we obtain the correct momentum distributions with this temperature 
playing exactly the same role as the "usual" temperature. 

For the case of a many-body system like a nucleus where a nucleon 
trajectory is chaotic, a typical eigenfunction will be written as above. Thus, 
as the system evolves and shape of the nucleus changes, the eigenfunction will 
acquire a geometric phase (see Sec. 8). It is this phase which will eventually 
find a relation with the absorptive part of generalised susceptibility jig] . 



4. Semiclassical linear response 



To study the effect of external perturbation on a many-body system, one 



can make usage of linear response theory [34, 35, pQ ] . Also, one can imagine 
a fictitious external field to study collective excitations as if they were gen- 
erated due to a hypothetic field p0|] . In this section, we obtain an expression 
for generalised susceptibility and the frequency-dependent response function 
(also called as the polarisation propagator J36|) in terms of periodic orbits 



employing the Gutzwiller trace formula 26]. We assume that the dynamics 
of a single particle is fully chaotic. As mentioned earlier, the important case 
of mixed dynamics will be taken up in Section 6. The assumption of chaos is 
consistent with some recent numerical inverstigations for the case of nuclei 



The system is described by a Hamiltonian H which is disturbed by a 
field, F ext (t). The total Hamiltonian is 

H T = H -QF ext (t) (7) 

where Q is an observable, an example could be magnetization in the context 
of spin systems, or, an electric dipole operator in photoabsorption experi- 
ments. The response function can be written as the imaginary part of the 
dynamical susceptibility, 

x"(M') = ^{{Q(t),Q(t')] 
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/ 



|V^W). (8) 



The angular brackets denote the expectation value and the square brackets 
denote the commutator. Setting an initial time to zero and final time to t, 
we can re- write 

x"(t) = ±([Q(0),Q(t)}) (9) 

where 

Q(t) = ek kt Qe-l k \ (10) 

(...) denotes an average over the initial state of the system which is, for 
instance, the thermal state wherein 

w- ^ •-**<...) an 

where the temperature, T = 

Denoting the many-body eigenstates by & n and one-body eigenstates by 
(j) n , and after reducing the description of the many-body system to one-body, 
the propagator is 

x»-gi^iwi ^^ + ": + 'o' - < i2 » 

where 

" FD < £ ") = Ti^WTi (13) 

denotes the Fermi-Dirac probability of occupation number with fi as the 
chemical potential. Eq. (12) is the Fourier transform of the two-time corre- 
lation function. In (12), the operator q is the one-body operator which when 
taken for each body and direct-summed gives us the operator, Q. 

Now that the susceptibility has been reduced to a 1-body expression, 
we can treat the 1-body system semiclassically in terms of periodic orbits 
in an effective Hartree-Fock (or some other mean-field) Hamiltonian, h. For 
subsequent convenience, let us re- write (12) as 

X ~"H = Y,\^a\m b )\\ FD {e a )-p FD {e a + %u)]8(nu + e a -e b ) 

a,b 

--E^rE ^^(e a mam)\ 2 s(hu +e a - £b ) 

n=l ' a,b 

2Ti , n! 
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Here we have introduced the time correlation function f n (t) and its Fourier 
transform, f n (u): 

fin FD 



/n(w) = jdte^Ut) 



Qn FD 

= 2nh^-^(e a )\{tl> a \q\<k)\ 2 S(nLJ + € a -e b ). (15) 

a,b 

In fact, 

f d n v FD 

Ut) = j de°-^(e)C(e,t) (16) 

with 

C(e, i) = tr <5(e - h)q(t)q(0). (17) 
Finally, the propagator is 



n=l 



The classical dynamics in the three-dimensional one-body effective po- 
tential of a fermionic system will, in general, be chaotic. Now we apply the 
Gutzwiller trace formula to obtain an expression for the propagator. For 
this, we need the semiclassical expression for the time correlation, C(e,t), 
which we consider in the form : 

C(e,t) = tr S(e-h)X 

= tr 5{e - h)A + itr S(e - h)B, (19) 

where X = X(t) = q(t)q(0) is one-body operator. This operator is non- 
hermitian but can be decomposed as X = A + iB in terms of two hermitian 
operators : 

j X + X* A X-jtf 

A=-5-, B = ^- . (20) 

We need to re-express matrix elements of A (and B) overthe eigenstates 
of h. These matrix elements can be obtained by employing first-order per- 
turbation theory for the perturbed Hamiltonian h(X) = h + XA. Assuming 
the eigenvalue problem for h(X) to be solved, the matrix elements of A may 
thus be obtained in terms of derivatives of eigenvalues of the perturbed 
Hamiltonian, e n (A), with respect to A. 
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Each term of the correlation function C(e,t) can be expressed as 

1 A 
tr 6(e-h)A = 3=tr 



vr e - h + i0+ 
= -Qftr -^logfe-Zt- Ai + iO+) 

7T OA 

On comparing with the identity, 

Id - - 

- -— 3ftr logfe - /» - AA + »0 + ) = tr <5(e - - AA) 
7r oe 



(21) 



we have 



- - ON 
tr < 5( e _^ = _|_( e;A ) 



f^A), (22) 



(23) 



where A^(e; A) is the cumulative density of levels for the parametrised Hamil- 
tonian h(X). 

For the case under consideration where the single-particle dynamics is 
chaotic in an effective mean-field whereupon we can assume that the periodic 
orbits are isolated and unstable, the cumulative density of levels is given by 
the well-known expression due to Gutzwiller |2q|, 



{2irh)f 

+ E yJ_^[I^A)-r|^] + 
p r =!™ |det(m;(A)-l)|3 

Here, irip(A) is the monodromy matrix governing the stability of the classical 
periodic orbits, p; v v denotes the Maslov index of the trajectory, and h\y is 
the Weyl symbol of the Hamiltonian operator, h(X). 

Owing to (24), we need the derivative of the action w.r.t. the parameter, 
A [071: 



dS p _ I dt dhw^ = _ , dtAw (25) 



d\ J d\ Jp 
where Aw is the Weyl symbol of A. Thus 

d* xd$ p 



tr Me - h)X = J ^LJLx w 6{e - h cl )\ + 0(T,-' +1 ) 

+ ^ E t ^ S '~ ^ / **w + 0(h°), (26) 
nn> „ __i I det(mnA) - l)| 2 J P 



p r=l 

where Xw is the Weyl symbol of X. 
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For our case where 



X = exp ( ^htj gexp f — —ht ) q, 



(27) 



the Weyl symbol can be written as 

AV(x,p) = (etV M Mi^ 



tiA 



{e -iht )we - 2 nA gw 



khA. 



exp(— Cdt) qw + Weyl-Wigner corrections 

= qw(t)qw{0). (28) 

The operator, C c i is the classical Liouvillian operator defined by {h, .}, 
and the operator A is defined by 



_ d d d_ d 

dp dx dx dp 



(29) 



and qw(t) = <7w[3>*(x, p)] with the Hamiltonian flow denoted by $*. 
Finally, the correlation function is 



C(e,t) 



tr 6(e - h)q(t)q(0) 
d^xd^p 



qw 



J_ v ~ cos [%S. 

p r=1 , , A ... p 



exp(-£ d i) qwS[e-h cl ] + 0{h' f+1 ) 
dTq w (T)qw{t + T) 



P r 2 V V\ 



f^i I det(m;(A) - I) 



2 JJ) 



+ (30) 
The propagator or frequency-dependent response function is given by 



dt 



d-> xd* p 



p FD (e) 



qw(*, p) 



\e y cu qw 



(x,p) + 



de 



p FD (e)- P FD (e+^) 



^g cos[§g p r\v p \ j 



p r- 



^|det(m;(A)-i; 



drq w (r)q w (t + Tpl) 



Since the rate of energy dissipation is directly related to x"i u )> the above 
expression gives us a semiclassical description of dissipation. Note that the 
shape (which means potential) of a system dictates the periodic orbits and 
the Liouvillian operator. 

5. Random- matrix linear response 

Energy level sequences of nuclei have local fluctuation properties which 
fit remarkably well with the predictions of random matrix theory [^] . So is 
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the case with other chaotic quantum systems | |32| [. An attempt to arrive at 
a connection between random matrix universality and chaos has been made 
recently [41]. In this section, we present random- matrix expression for linear 
response for Fermi systems. Thus the Fermi system is one in an ensemble, 
where the Hamiltonian matrix is modelled in terms of a random matrix. In 
the context of absorption by small metallic particles, random matrix theory 
was used by Gorkov and Eliashberg [42]. 

Given a level sequence, {xi}, the n-level correlation function is defined 
as @ 

Nl 

R n (x 1 ,x 2 , ...,x n ) 



(N — n)\ 

dx n+1 ...dx N P N (x 1 ,x 2 , ...,x N ) (32) 

which is the probability density of finding a level (regardless of labelling) 
around each of the points x±, x 2 , x n , the positions of remaining levels 
being unobserved. P{x\, ...,xjv) denotes the joint probability distribution 
function for the levels x\, xjy. The level density is given by the one-point 
function, Ri(x). The n-level cluster function is defined as 

m 

T n ( Xl ,x 2 , ...,x n ) = ^(-l) n ~ m (m-l)! J] R Gj (x k , with x k e Gj),(33) 

G 3=1 

where G stands for any division of the indices (l,2,...,n) into m subgroups 
(Gi,G 2 ,...,G m ). For instance, 

T l (x) = R 1 (x), 

T 2 (xi,x 2 ) = -R 2 (xi,x 2 ) + R 1 (x 1 )Ri(x 2 ), (34) 

and so on. Since we are only interested in the local fluctuations, the energy 
levels Xi and Xj are such that N\xi — Xj\ ~ 0(1), the so-called scaling 
limit. In the scaling limit, the two-level cluster function, T 2 (x±, x 2 ), becomes 
Y 2 (\ Xl -x 2 \). 

Time correlations will now be averaged over the ensemble of random 
matrices which entail an ensemble-averaged time correlation function. Each 
term is multiplied by the probability that two levels will be in the intervals 
dei and de 2 , 

l2(|ei - e 2 |)^^ = y 2 (|ei - e 2 \)n 2 T F de x de 2 (35) 
a a 

where a is the mean level spacing and utf is the Thomas- Fermi density of 
levels. We get for the response function : 

2 ^ (M m f°° f°° j j d m p FD (e) ... 
X» = n 2 TF £ i-f- / / de x de 2 -e 2 + e 1 ) 

m=l J —oo J— oo ft 

l2(|ei-e 2 |)|(0i|^ 2 )| 2 . (36) 



12 



For the invariant ensembles, i.e., the orthogonal, the unitary, or the sym- 
plectic ensembles, — £ 2|) embodies universal results. Thus the random 

matrix expression gives the generic form. It must be borne in mind that 
the local fluctuations and universalities have been qualitatively (and some- 
times quantitatively) understood using semiclassical methods. The random 
matrix result provides us a trend, a qualitative overall behaviour; it is the 
semiclassical theory containing microscopic information that can eventually 
give us system-specific results. 

The matrix element appearing in (36) can be related to the classical 
autocorrelation function of q c i (the classical counterpart of q) if we assume 
that the sum of mean values of any quantity over all states with definite 
energy roughly equals the sum of the classical mean values of this quantity 



as the particle moves over all the trajectories with given energy [50]. 
Thus, 



m=l 



2tt 



ml 



m^FDi 



de\de\ 



d m p 



de n 



/oo 
dre^-^ T C gq (r) 
-oo 



where the correlation function is 
1 rT 



C qq (r) 



lim 

T^oo T 







dtq(t)q{t + T). 



-8(Tuu - e 2 + ei) 
(37) 

(38) 



Changing the variables to 



£2 - ei 
£2 + ei 



and after some manipulations, we get a compact result : 

x"(uj) = nTF^YY 2 (huj)S q q(Lj). 

In this expression, S qq {uj) is the spectral density defined by 



SM = - 



dTe^ T C qq (r). 



(39) 



(40) 



(41) 



For different Gaussian ensembles [4C] and non-Gaussian, invariant ensembles 
[43 1, the two-point function is 



y 2 (M 



sin(-7r7iu;) 



n 2 



sin(-7r7iu;) 



-i 2 



sm(2irhuj) 



+ 



/ sin(7rt) 



dt 



d ( sin(7ri) 



d(fiuj) 



sin(27r?ia;) 
27rhtu 



d(hu) V 7ri 
hu rsin(27rt) 
,, _ 2-Kt 



dm 
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corresponding to unitary, orthogonal, and symplectic ensembles respectively. 

For the case when a random matrix, h, is perturbed by another random 
matrix, q, independent from h but invariant with respect to the same symme- 
try, S qq (uj) is explicitly known |5j|. The spectral density becomes asymmetric 
at non-zero temperatures and shows a universal dip at zero frequency. The 
behaviour of the spectral density near zero frequency is universal, depending 
on the symmetry parameter. 



6. Response of mixed systems 



The most common situation is when the dynamics is neither fully in- 
tegrable nor fully chaotic, but mixed. Looking at the phase space of such 
systems, we observe stable islands embedded in stochastic sea. We present 
an expression for the response function employing the recent semiclassical 



result for mixed systems [16|. 



The semiclassical expression for the density of energy levels is given by 
the trace of the Green function in position representation, and it results in a 
sum over periodic orbits. The integration in the method of stationary phase 
goes from one point x to itself collecting contributions of periodic orbits. 
Here the symmetry becomes important. If there is a continuous symmetry, 
then there is a family of periodic orbits with given energy. Only the inte- 
gration transverse to this family can be done by the method of stationary 



phase, the parallel integration remains. Following [16], let us assume that 
the symmetry group is G. In case G acts on the phase space in a way that 
infinitesimal generators of the group, and the vector field of the Hamilto- 
nian itself, are all linearly independent, the subgroup of G which leaves the 
periodic orbits invariant should be discrete. In this case, the degeneracy of 
periodic orbits , k is just the dimension of the group. Also, if 70 (i) is a refer- 
ence orbit, other orbits can be parametrised by group elements g according 
to 

lg {t) = g lQ {t). (43) 

The points on the family are parametrised by (g, t) and the (fc+l)-dimensional 
measure is dtdfx(g) where dfi(g) is the Haar measure |Q. The trace formula 



for the case of exact symmetry becomes [45], 



1 1 » — , f 1 iS(E) icTTT 



po,r 



where K is an invariant of the family and is determined by linearisation of 
dynamics about a typical orbit. 
For nontrivial symmetry, 

K = QdetW det(M - I). (45) 
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Here Q is the jacobian, independent of dynamics from the choice of basis 
for the Lie algebra, W is the symplectic matrix from the linearisation in full 
phase space of the dynamics around the orbit, and tilde on M represents 
the fact that symmetry has been extracted. In (44), S(E) is the action of a 
typical orbit in the family, and, a = fi — 5 where \x is the Maslov index and 
5 is the number of positive eigenvalues of W. 

Now we perturb a system (with Hamiltonian Hq) with continuous group 
G so that the resulting Hamiltonian, H has no symmetry at all, 

H = H + XH 1 . (46) 

For this case, the trace formula becomes II? 
1 1 



where 



e— -(e^ F{9 ' E) ) geG , (47) 

f AC 

F(g,E) = - dtH 1 = —{x ]] ,E;\). (48) 
J lg {t) A 

The angular bracket represents an average over the orbit label g. S$ is the 
action of a trajectory of unperturbed Hamiltonian and AS is a correction 
of order A. To provide an expression for two-time correlation function, we 
first note, progressing in the same way as in Sec. 4, 

5Jt*w»» 1 A 5 e " ^i»(x E ) (49) 

where 

M(j,E^=(A F ^% eG , (50) 

and X is defined in Sec. 4. In general, let X be A(0)i?(i). For the case of 
mixed dynamics, we obtain the following expression where the system is in 
thermal equilibrium at temperature, T: 

C AB (t,\,T) = ^ r J'^2 e -6»..A w (x,p)e- c -'B w ( X ,p) 



iSg(E) ia7T d / 
en~, _/exp 



— / dr A(r)B(r + t') 



'(51) 

J / geG 
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Once the group G is identified, the calculation of <> g eG 1S n °t difficult. 
For the case of axial symmetry, it has been demonstrated in [jl(J. In the 
case of nuclear physics, the connection between deformed nuclei and har- 
monic oscillators where the frequency ratios can be related to deformation 
parameters, a much desirable simplification is expected to occur. 

Note that the decay of correlations will be governed by the eigenvalues 
of the classical Liouvillian operator. For the systems which are mixing, the 
decay of correlations is exponential. For mixed systems, this decay will be 
more complicated and system-specific. The important point is that dynam- 
ical aspects have been incorporated explicitly. We believe that it will lead 
us to an understanding of many important aspects of many-body systems, 
pre-fission neutron multiplicity, damping of giant resonances being some of 
the examples. 



7. Irreversibility of energy dissipation 



It has been argued in the past that the irreversibility of energy dissi- 
pation occurs when the single particle dynamics is chaotic ||. Quantum 
mechanically, the difficulty is in showing that the energy distribution, rj(E), 
defined below follows an equation which implies irreversibility. Most often, in 
literature, this equation is taken as the Smoluchowski equation. Within the 
assumptions of classical linear response, this is the correct equation, as shown 
by Jarzynski [46|. However, our system is quantum mechanical though we 
have simplified our considerations somewhat by assuming that our system 
is adiabatic. In the two subsections below, we derive evolution equations 
for the energy distributions and show that these are, indeed, irreversible 
in time. We assume that the single-particle dynamics is non-integrable. It 
has been shown for chaotic quantum systems and also for non-integrable 
systems with zero Kolmogorov-Sinai entropy that there is level repulsion 
[32, 33]. That is, the probability of clustering of two adjacent levels is zero. 
In the quantum case, the equation contains a term in addition to the Smolu- 
chowski equation. The subsection on semiclassical adiabatic systems shows 
that the quantum diffusion equation, obtained by Jain |47|1 , reduces to the 
Smoluchowski equation as h — > 0. 



7.1 Quantum adiabatic systems 



We consider a quantum system evolving adiabatically in time. Since the 
evolution is adiabatic, there is an instantaneous basis. Thus, the energy levels 
evolve in time as the system "deforms". While in the same symmetry class 
(we assume so), the levels keep evolving without crossing. The instances 
when the levels come closer than the mean level spacing, there is non-zero 
probability of transitions to take place. For instance, very small spacing be- 
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tween levels then leads to an increased probability of non-adiabatic Landau- 
Zener transitions |3l]] which eventually, in long time, modify the energy 
distribution of the system. The basic idea behind using the Landau-Zener 
transitions goes back in the literature of nuclear physics where it was 
used to argue for damping of collective modes. 
Let us consider the Hamiltonian, 

H{t) = H + etV (52) 

where Hq,V are linear operators. We assume that the evolution in time is 
adiabatic which corresponds to the smallness of e. 

At any instant, the system admits an eigenvalue spectrum given by the 
eigenvalue problem for the "frozen" Hamiltonian, 

H(et)\n(et)) = E n {et)\n{et)). (53) 

From time t = 0,the levels evolve in time, the resulting density operator, 
p satisfies 

ih^ = [H,p]. (54) 
Our objective is to derive an equation for the energy distribution, 

r](E) = / dE'tr{5(E' - H)p}. (55) 

When a system is perturbed adiabatically, there is, of course, a clear sep- 
aration of time-scales. To incorporate these scales in the problem, we employ 
the multiple time-scale method for treating the partial differential equation 
(54). Accordingly, denoting the set of instantaneous states by {\n(et))}, we 
can write an expansion for the density operator, 

p({\n(et))}, t) = Po ({\n(et))}, et) + e Pl ({\n(et))}, t, et) + ... (56) 

with the initial conditions, 

po = p o(H(et)), pi = p 2 = P3 = ••• = 0. (57) 

Substituting (56) in (54), we get a system of equations separated by different 
orders of e : 

\po,H(et)] = 0, (58) 
+ = -^fe, J = 1,2,... (59) 

If there are no other constants of the motion than H(et) on the fast scale, 
or under the Thomas- Fermi approximation, by (58), 

p ({\n(et)}},et) = p' (H(et),et) (60) 
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where the arbitrariness of p' is removed by insisting that p remains valid for 
times 0(e~ 1 ) by removing secularities in (59) with j = 1. 
To remove the secularities in 



ih— + [ Pl ,H(et)] = -zh m 



(61) 



in keeping with the method, we operate by an arbitrary operator-valued 
function, g(H) and perform the trace of the resulting equation over the 
frozen basis, 



/ dpi 



n 



ih 



d(et) 



n 



(62) 



We should finally regard those results which hold for any g(H). 

For pQ to be valid for times for 0(e~ r ), the right hand side (RHS) of (62) 
should be set to zero, which leads to 



dp' (H) dH dp' (H) 
dH d(et) d(et) 

Let us re-write the density of levels as 
£(£7,ei) = ^(n|5(E-^)|n) 

n 

d 



n 



0. 



(63) 



dE 



J2{n\&(E - H)\n) 



dnjE, et) 

dE 



(64) 



where is the cumulative density of levels. In the following, it will often 
be useful to convert the traces into energy- integrals. For that, we define an 
energy average over the "frozen" Hamiltonian by 



^ n 

With this definition, it is easy to verify that 
£>|...|n> = f dEZ(...) E , et . 

n J 

Now, (63) becomes 



(65) 



(66) 



S 

Calling 



dp'o I dH 



dE \ d(et) 
' dH \ 



E,et 



dp'o 
d(et) 



0. 



(67) 



u(E,et), 



(68) 
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and using (66), we obtain the quantum adiabatic theorem: 



d(et) dE 



(69) 



In the sequel, the notation d x stands for a partial derivative with respect 
to x. 

This is a continuity equation for the level density. If the energy levels are 
thought of as fictitious particles, the level density is like the particle density. 
Indeed, in random matrix theory, one can imagine that the levels are like 
particles with Coulomb interaction in two dimensions [40, 55]. 

Therefore, (67) reduces to 



Then, for po> we have 

wj mu) = ~dE ( 1 ] \m 



(70) 



(71) 



With this equation and the initial condition (57) , we have completely deter- 
mined pQ. 

Similarly, we can proceed to determine p\ [47|. 

Let us define the average two-time correlation function, 

C et (s, E) = ({d et H({\n)}, et) - u}{d et H{{\N}}, et) - u}) E ^ t 

= ±J2( n \ S ( E -H){detH({\n)},et)-u} 



(72) 



{d €t H({\N)},et)-u} \n 

and a coarse grain over the energy spectrum, replacing thereby 
g{E n ) by g(E), and 



n 



9 2 P' 



dE 2 



n ) by 



dE 2 ' 



Also, for times of 0(-), 

r° l r°° l 

J t dsC{s) = 2j.J sC{s) = 2° 2 



(73) 



(74) 



Finally, the condition that removes secularities to 0(e 2 ) is 




dp/A i^dp' dG 2 



2 dE dE 



0. (75) 
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The energy distribution follows the following equation, 



dr] d . . 2 d 

dt=- e dE {ur])+e dE 



4^(1). pi 



the quantum diffusion equation. Thus, the diffusion in quantum systems has 
to be qualitatively and quantitatively different as the diffusion coefficient will 
be different from the one we have in the Smoluchowski equation. 

There are different important time-scales in the system. First of all, the 
time scale associated with the decay of correlation function, 

/+oo 
C(s)ds. (77) 
-oo 

If the quantum system considered is modelled as a random matrix of di- 
mension N |53| with large N, we know that correlation function will decay 
very rapidly. Thus, t c can be very small if the quantum systems possess the 
following properties : (a) the number of eigenvalues is very large, and the 
energy spectrum is complex, and, (b) the corresponding classical system is 
chaotic. Chaos in the underlying classical system plays a fundamental role 
in the decay of correlation functions. It was shown by Gaspard and Jain 
[54 1 that the quantum time-dependent correlations in a Fermionic system 



are dominated by the classical correlation function. The decay of the cor- 
relation function is shown in this work to be governed by the eigenvalues 
of the Liouvillian operator. Thus, t c is related to the Liapunov exponents 
and other detailed features of chaos. In classical ergodic adiabatic systems, 
the time t (fast scale) is much larger than t c , thus the third term of (76) is 
zero. However, in quantal systems, we have the quantum mechanical scale, 
t q = h/W (<r being the mean level spacing) which is why the third term at 
0(e 2 ) is explicitly present. If t q <C t c <C t, the quantum effects will dominate, 
and all the terms in (76) will be important. If t c <C t q <C t, then the sys- 
tem will behave classically initially and eventually, quantum phenomena will 
become important; so initial evolution will be Smoluchowski-like and then 
non-Smoluchowski regime sets in. If, however, t c <C t <C t q , then the evolu- 
tion will be according to the classical equation. Notice, as Ti becomes small 
and the system is classical, t q will become vanishingly small. In the limit, the 
correlations will become largely independent of energy and the third term 
of (76) will drop. We now show that the equation in the semiclassical limit 
is indeed the classical Smoluchowski equation. 



7.2 Semiclassical adiabatic systems 



The combination of semiclassical and adiabatic approximations presents 
a singular situation. We present here a systematic treatment where we per- 
form an ^-expansion in the Wigner-Weyl basis on top of the multiple-time 
scale expansion. Proceeding in the same way as for the quantum adiabatic 
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case, we will derive an equation for the energy distribution. The density of 
levels is defined as 



E(E,et)=J2(n\S(E-H)\n). 



(78) 



In the Wigner-Weyl basis, it can be written through the Wigner transform 
ol5{E-H) H], 



8(E - H) w 



5{E-H) + Ti 2 \ -5"(E-H) 

2 



1 



6"'(E-H) — 



l d 2 V 

2 



Here, the Hamiltonian is assumed to be of the form ^- + V(x), and the 
subscript W refers to the Wigner transform. The Wigner transform of X is 
consequently written as 



dpdx5(E - H) 



J dpdx5(E -H)-— J dpdx5"(E - H)-^ 



(80) 



The primes on the Dirac delta distribution denote the weak derivatives or 
the distributional [49]. The energy average of an observable, Q, is 



< Q >E,e 



AY 



dpdx5(E - H) W Q. 



(81) 



The Wigner transform of the density operator used in the earlier sub-section 
is the Wigner distribution, fyy. The evolution equation for the Wigner dis- 
tribution is obtained by taking the Wigner transform of the von Neumann 
equation : 

dfw 



at 



4 sin Q A ) /w 

// ( A - ^A 3 + ...)f w . 



(82) 



Upto 0(h 2 ), for the Hamiltonian of the form specified above, the evolution 
equation for fw is 



dV df w df w ft 2 d 3 V d 3 f } 



dt <9x dp ^ dx 24 <9x 3 dp 3 
We now expand fw in a power series of Ti : 

fw = f {0) +h 2 fW+h*f^ + ... 



(83) 



(84) 
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As a result, we obtain the following hierarchy of equations : 
0(h°) : d t fM = {H,fW}, 



24 

0(?i 2n ) : 9 4 /( 2n ) = {F,/ (2n) } - ^dlV.dlf^ n -^ . (85) 

{., .} denotes the Poisson bracket in the mock phase space variables. To solve 
these equations, we shall exploit the two time-scales, t and et, and expand 
/(°) ((x, p) = z G mock phase space |^8| ) as 

/(°) = / (0) (z,et) + e/ 1 (0) (z, e M) +e 2 fi°\z,et,t) + ... (86) 

with the initial condition, 

/,S 0) (z,0) = /oo, ff ] = V i > 1 at t = 0. (87) 

To various orders in e, we have a hierarchy of equations nested in the fi- 
hierarchy: 

O(e ) : {/ o (0) ,iT} = 0, 

(0) a/ (o) i 



0(0 : ^T+{fi°\H} 



dt lJn ' J a(et) 

Due to the O(e ) equation, 

f(°\z,et) = f (H,et) (89) 
with the initial condition, 

ME,0) = ME). (90) 
/o is completely specified by (90) and 

(> (S/ ) + A (nS/o) = (91) 



d(ei) v JUJ OE 

removing secularities at O(-) in time. In (91), we have employed the follow 
ing notations : 

d 



Z(E,et) = J 'dz5(E-H) = —n{E,et), 
l -Jdz5{E-H)..., 



< ••• >E,et ~- ^ 



/dH\ 
\d(et)/ 



(92) 
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Eqs. 



, (91), and (92) imply 
(°) / QH 



The formal solution of 



d(et) 
for n = 1 is 



it(i7, ei) 



* <9f (0) 



o 



(Z,et) + h(H,et), 



(93) 



(94) 



where Z = Zi(z,t,t' ,et) is a point in mock phase space reached by starting 
at z at time t, then evolving backward in time to t' under H. 

We can determine f[°^ by removing secularities at 0(e~ 2 ) of time. Before 
doing so, we need some definitions which will be encountered in the sequel. 
The two-time correlation function is 



c(.) 



oh 



u Of*(s. 



f dH 
\d(et) 



(95) 



where the operator acts to its right, evolving a point z for a time s under 
H. C(s) satisfies the same conditions as in the previous sub-section. The 
condition avoiding secularities at 0(e 2 ) is 



_|_ (E/l)+ > E/l) _^( EG2 §§ 

Thus, valid to 0(e), the distribution function is 
/ {0) (z,et) 



0. 



(96) 



fo{H,€t) + efi{H,eb) 



o 



\d(et) 



[Z,et) -u 



(97) 



where /q and f± satisfy (91) and (96) respectively. 

Now we expand the quantum correction, f( 2 \ as /(°) with fi and e as 
independent parameters : 



= /f(z,et) + ef?\z,et,t) + e 2 (z,et,t) + .. 
To orders in e, we have 



(98) 



0(1) : {H,f§ i >} = ^%.%f l 



df 



(2) 



ot 



The solution of (99) (for / (2) ) is 



-a 3 a 3 f. (0) 



/o (2) (z,e*) 



2 AO) 



d(et) 24 



±(d x V) 2 + ±(p.d x fV(x) 



(99) 



a|/ (0) }.(ioo) 
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Upto 0(e) and 0(h 2 ), the semiclassical density is 

f w (z,et,t) = fo{h,et) + efi(h,et) - e— (h,et) 



dt' 



o 



dh 



(Z, et) — u(h, et) 



4/f \. (ioi) 



i(d x U) 2 + i(p.d x )V(x) 

Now we consider (99). Multiplying by a test function, g(h) in the equation 
for % = 1 and integrating over the phase space gives us 



d_ 

dt 



J dmg(h)f® = -J *9{h)^- -±J dk^)fl£V.^/{ 0) . (102) 



Since the RHS is independent of t, the LHS grows secularly. We set the RHS 
to zero, and after some manipulations, we get 







(S/ (2) ) + j^VZfn + J dEZ ^ = 0. 







9(6t) 



1 



r(0)\ 



(103) 



With (91), (96), and (103), we obtain the semiclassical diffusion equation, 

d ( 7/ 



<9// 

at 



9 , x 1 2 s 
e— (u77) + -e 



as 



2 9^ 



G 2 S 



5S VS 

3 f (0) 



(104) 



E,et 



(2) 

We have not succeeded to include f\ in this work. 

Note that (104) reduces to the Smoluchowski equation as h — > 0. 

8. Geometric phase and dissipation 



Adiabatic approximation leads to linear response theory on one hand 
where dynamical susceptibility (or polarization propagator) is central, and 
geometric phases on the other. When a particle (e.g., a nucleon) moves inside 
an enclosure whose boundary is adiabatically vibrating in time, the wave- 
function can acquire a geometric phase over a cycle of vibration. In this sec- 
tion, following flql , we show that the geometric phase thus acquired is related 
to the imaginary part of susceptibility. It is well-known that the imaginary 
part of susceptibility is a quantifier of dissipation, geometric phase plays an 
important role in understanding of damping of collective excitations. 

To understand how this relation comes about, for the sake of brevity, 
we restrict ourselves to the case of cyclic evolution of deformations of the 
many-body system. Let us denote the Hamiltonian parametrised by some 
parameters (e.g., deformation parameters), R, by H(R) which describes a 



single particle in an effective mean- field. It is well-known [59| that when 
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the parameters evolve along a cyclic path C the instantaneous eigenfunction 
of the system | n(R)) corresponding to the eigenvalue, E n (R), acquires a 
geometric phase given by 



7„(C) = £i(n(R) | V R n(R)).dR 

V n .dS, (105) 



/c 

1 

~h 

where S is the surface enclosed by C in the parameter space, and V n is 
the "field strength" (adiabatic curvature) given by a familiar expression 
involving a wedge product : 

V n = -ih ^ — — — . (106) 



A suitable form of V n for the sequel is [61] 



V n = — lim / dte~ et t(n\ [(V R H) t , A(V R ff)] \n) (107) 

where (V R iI)i denotes the Heisenberg-evolved operator. Note that, the state 
|n(R)) appearing in (105) corresponds to a single-particle eigenket in an 
effective mean-field. This state is clearly related to the original many-body 
Fermi system for which the imaginary part of the dynamical susceptibility 

is mil 

X"(t) = ^<*oM(f),£(0)]|*o> = / ^e--'x>) (108) 



where |3>o) is the pure ground state of the many-particle system with Fermi 
energy, Ep. If one-body operators, H, A, and B, are used to construct many- 
body operators by a direct sum so as to get 7i, A, and B, respectively, and 
H^) = £i\^i), then we |3 have 



($ \A\$i) = <m(R)|A|n(R)). (109) 

On reducing the many-body system at T = OK, (where the Fermi-Dirac 
distribution is a Heaviside step function), to one-body system, we can ex- 
press [54] 



= -|/ dte^ tr 5{E F - H)[A(t), B(0)}. (110) 

This can be written as in Sec. 4 semiclassically. 

The label n in (105) corresponds to single-particle states and is related 
to |<J>o) because the many-body matrix elements can be written in terms of 
one-body matrix elements for the case when all the constituents are taken 
as non-interacting. In many-body physics, this gives the zero-order response 
whereupon the interaction can be included in a Vlasov description in an 
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iterative way |)3|]. For relating the response function to the geometric phase, 
the operators A and B in our discussion are to be identified with VjW and 
V Y H for R= (X,Y,Z). 

The matrix element in (105) can be written as a many-body matrix 
element using (109) by composing A and B so that we get the operator, 
C(t) = [A(t),B(0)] - [B(t),A(0)], which is related to a difference \" AB {t) - 
Xba(^) = Xc(f)- Thus, we can re- write V n as 



V n = rrlim/ dte~ et t (* |C|* ) 







dt t X 'c(t) 



■jj=0- 



duj 

We now arrive at our relation for the case of cyclic evolution : 



7„(C)= jf dS. 



(111) 



(112) 



Since Xc( u; ) ^ s related to energy dissipation, this relation connects geometric 
phase to dissipation. In other words, we see that dissipation in finite systems 
is purely quantum mechanical as it is related to geometric phase. That this 
forms a fundamental basis for the Bohr's liquid drop model is discussed 
elsewhere |l5[| . 

Similarly, we can find a relation for the case of non-cyclic evolutions. 

It is known that friction or viscosity in a quantal system can appear 
in a thermodynamic limit. Since we have a many-body system, seemingly 
contradictory conditions of finiteness of the system and "continuity" of the 
energy spectrum are met with. The finite size explicitly manifests itself in 
terms of a sum over periodic orbits. 

It is important to note that the wavefunction that we have considered 
above is the one which satisfies the requirements in Section 3. 

It is interesting to note that the viscosity of quantum Hall fluid in two 
dimensions at zero temperature is related to adiabatic curvature [S2|. Also, 



for covalent dielectrics, the polarisation is related to geometric phase. These 
two examples and our general treatment leads us to think of a deeper unified 
connection. 



9. Fission viscosity tensor 



Since long, linear response theory is used to study the dynamics of fis- 
sion process p5| in nuclei. The collective coordinates are the deformation 
parameters and pairing gap. The quantity related to fission viscosity tensor 
is the first moment of the time-dependent response function. In this pa- 
per, we have presented semiclassical and random matrix expressions for the 
response function. In this Section, we apply this knowledge to obtain expres- 
sion for viscosity tensor in terms of periodic orbits of the underlying classical 
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system. It will also turn out that the classical Liouvillian operator plays a 
fundamental role. This is very interesting as the eigenvalues of the Liouvil- 
lian dictate the decay of correlations. Thus, the existence of dissipation is 
directly related to the non-integrability of the classical system. 

If we denote by {q u } the set of N collective coordinates, then the friction 
viscosity tensor is given by 



C UfJ ,(t,T). The system is at a temperature, T, related to the excitation en- 
ergy of the intrinsic system. It was brought into the calculations when the 
intrinsic degrees of freedom were averaged over the canonical ensemble. Once 
again, we recall that the canonical statistical mechanics follows from the as- 
sumption of microscopic chaos in single-particle motion. 

The semiclassical form of the time-dependent response function (or the 
time correlation function) gives the following expresssion for the viscosity 
tensor when we assume that the single-particle-dynamics is fully chaotic (in 
case the dynamics is mixed, expression can be immediately written using 



This is a semiclassical expression based on incorporation of the Gutzwiller 
trace formula in the linear response theory. Applications to the practical sit- 
uations in nuclear fission will be of a great interest. 

A random matrix expression for the viscosity tensor can also be readily 
written employing the results of Section 5. The coefficient of friction in the 
zero-frequency limit is 



where is given by (40). Notice that friction is physically the adia- 

batic curvature, from (112). This has been also found earlier for viscosity of 




(51)) : 





(115) 



quantum Hall systems |62| . 



10. Summary 
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A lot of work has been done on the nature of nuclear dissipation (see, 
e.g., HH). It has been stated, most clearly by Weidenmuller ||, that chaos 
in the nucleonic motion will be related to dissipation. However, before trying 
to explain the experimental data, we need a systematic theory. In this work, 
we have taken a step in realising this connection. 

On the basis of the connection established between quantum chaos and 
statistical mechanics, we have developed semiclassical and random matrix 
theory of response to external perturbations of a many-body Fermi sys- 
tem. For the cases where the dynamics of a single particle is chaotic or not 
fully chaotic but mixed, we have presented the results for response func- 
tions. Also, we have shown that the collective excitations damp irreversibly 
and this irreversibility of damping is shown to be quantum mechanical. The 
damping of collective motion is reflected in the re-distribution of energy. 
Energy re-distribution has been shown to evolve irreversibly. The reason for 
this irreversibility is the complexity of the energy spectrum which might, in 
turn, be connected with classical non-integrability (chaos or mixed dynam- 
ics). The complexity is such that, in the scaling limit of the energy spectrum, 
the system behaves as if it is a realisation of an ensemble of random ma- 
trices. Quantum chaos or complexity of spectrum brings with itself a kind 
of almost periodic evolution so that the system takes much longer to show 
any recurrence. The recurrence is a must for a system with a discrete spec- 
trum. However, collective excitations decay much before this recurrence can 
occur. To present a simple instructive example, we consider the set of prime 
numbers and look at 

V{t)= ]T esxp[ipt], (116) 

primes, p 

which can be thought of as the survival probability (a special case of two- 
time correlation function) of a fictitious system. This sum decays quickly as 
the number of primes increase. Clearly then, for a many-body sytem like a 
nucleus where the spectrum is quite complex, the time correlations decay 
very quickly. This will be decided by the spectral properties of the Liouvillian 



[ 28 1 . The coarse-graining we had done in Sec. 7 was a representative of the 
complexity we have discussed here. 

Another point of importance is the arbitrariness of the initial density 
operator. We have obtained quantum and semiclassical generalisations of 
the Smoluchowski equation. Damping is described in terms of the imaginary 
part of the response function. This is shown to be related to geometric phase 
acquired by a single-particle wavefunction as the many-body system evolves. 

The application of the general theory based on existence of chaos or 
mixed chaos in the many-body system is being applied to fission of nuclei 
and metallic clusters. We have presented an explicit expression for the fission 
viscosity tensor in Sec. 9 which can be used to analyse the fission data for 
relatively heavier nuclei. Clearly, the quantum diffusion equations obtained 
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here are directly applicable to fission processes [|56], [57j where the Smolu- 
chowski equation has been used until now [58]. The difference new term will 
make in the assignment of transport coefficients is of interest. In case of nu- 
clei, in the past [67], there has been a criticism in modelling the deforming 
mean-field by the form, QF ext (t). However, this is basically a matter of con- 
venience. However, while making realistic calculations, one may keep this in 
mind. The advantage one has here is that one gets the response function. 

With the semiclassical wall formula obtained here, we are now trying to 
implement it by finding periodic orbits of the underlying classical system. 
Since the trace formulae can be written in terms of zeta functions, the nu- 
merical implementation will be efficient. Hopefully, not too many periodic 
orbits will be needed. In this regard, since the deformed nuclei can be mod- 
elled in terms of harmonic oscillator potentials, and we have the exact trace 
formulae for them [j6"6f |, it will be interesting to use our expressions to find 
response functions. 

We hope to study nonlinear response in a semiclassical and random- 
matrix setting in future. This is important for understanding phenomena at 
an ultra-fast time scale. 

Hence, we have shown that there is an emergence of a new understanding 
of damping of collective excitations, on quantum transport coefficients in 
finite Fermi systems, and in the general theory of response of nuclei, metallic 
clusters, and quantum dots. 
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